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1 Introduction 

Nonlinear dispersive waves are wave phenomena resulting from the interacting effects of nonlinearity 
and dispersion. Dispersion refers to the property that waves of different wavelengths travel at dif¬ 
ferent velocities. This property may, for example, be due to the material properties of the medium, 
e.g. chromatic dispersion [lj, or to the geometric arrangement of material constituents, e.g. a 
periodic medium with Floquet-Bloch band dispersion [3]. Nonlinearity distorts the shape of a local¬ 
ized structure or creates amplitude dependent non-uniformities in phase. This may concentrate or 
localize energy in a region of space (attractive / focusing nonlinearity) or tend to expel energy from 
compact sets (repulsive / defocusing nonlinearity). In electromagnetics, nonlinearities may arise 
due to the intensity dependence of dielectric parameters (e.g. Kerr effect [I])- Physical phenomena 
in which the effects of both dispersion and nonlinearity play a role are ubiquitous. Some examples 
are: (a) long waves of small amplitude at a water-air interface [ 112] . (b) a nearly mono-chromatic 
laser beam propagating through air, glass or water [36], (c) light-pulses propagating through optical 
fiber waveguides [TJ and (d) the macroscopic dynamics of weakly correlated quantum particles in 
a Bose-Einstein condensate; see, for example, [8611301149| . Interest in nonlinear dispersive waves 
and their interaction with nonhomogeneous media ranges from Fundamental to Applied Science 
with great promise in engineering / technological applications due to major advances in materials 
science and micro- and nano-structure fabrication techniques; see, for example, m- 

An important feature of many nonlinear dispersive systems is the prevalence of coherent struc¬ 
tures. They range from phenomena we can only passively observe in nature, such as large-scale 
localized meteorological events, to those we can control, such as the flow of energy in computer 
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chips and other nano-patterned media. In general, a coherent structure is one which has a spatially 
localized core, e.g. a solitary wave concentrated in a compact region of space or a front-like state 
which transitions between different asymptotic states across a spatially compact transition region. 

We focus on a key model in the mathematical theory of nonlinear dispersive waves, the nonlinear 
Schrodinger / Gross Pitaevskii (NLS/GP) equation: 


id t <$> = -A$ + I/(x)$ + 5|^>| 2 ^ • (1.1) 

In particular, we consider its nonlinear bound states. Here, <h(x, t) is a complex-valued function for 
and t 6 M. V(x) is a real-valued “potential”, and g is a real parameter. The following table 
lists the physical significance of <3?, V and g in the settings of nonlinear optics and macroscopic 
quantum physics: 



Nonlinear Optics 

Macroscopic Quantum Physics 

<3?(x, t ) 
V{x) 

9 

Electric field complex envelope 
Refractive index 
Kerr nonlinearity coefficient, g < 0 

Macroscopic quantum wave function 

Magnetic trap 

Microscopic 2-body scattering length (g > 0 or g < 0) 


1.1 Outline 

Basic mathematical properties of NLS and NLS/GP are discussed in section[2j In section[3]nonhnear 
bound states and aspects of their stability theory are discussed from variational and bifurcation 
perspectives. Examples are also presented and the particular cases of V (x) given by a single-well, 
double-well potential and periodic potential are discussed in detail. In section[4]we discuss particle¬ 
like dynamics of solitons interacting with a potential over a large, but finite, time interval. In 
sections [5] and |6] we consider the very long time behavior of solutions to NLS / GP. In particular, the 
focus is placed on the important resonant radiation damping mechanism that drives the relaxation 
of the system to a nonlinear ground state and underlies the phenomena of Ground State Selection 


and Energy Equipartition. Sections 5.1 5.2 discuss linear and nonlinear toy minimal models , which 
illustrate these mechanisms. Regarding the overall style of this article, we seek to emphasize the key 
ideas, and therefore do not present all detailed technical hypotheses, leaving that to the references. 


Acknowledgements: The author would like to thank J. Marzuola and E. Shlizerman for stim¬ 
ulating discussions and helpful comments on this article. He also wishes to thank Bjorn Sandstede 
and the referees for their careful reading and suggested improvements. 
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2 NLS and NLS/GP 

In this section we review basic properties of NLS/GP: 

idt$ = - A4> + P(x)4> + c/|4>| 2ct 4* . 


( 2 . 1 ) 


Here, g is a real parameter and we assume that V is a bounded and smooth real-valued potential. 


The equation (2.1) is often called the nonlinear Schrodinger / Gross-Pitaevskii (NLS/GP). 


We consider solutions to the initial value problem for (2.1) with initial conditions 


4>(x,0) = $ 0 (®) e H\R d ). 


( 2 . 2 ) 


The basic well-posedness theorem states that for d = 1,2 and all 0 < <r < oo and for d > 3 


and 0 < a < 2 /{d — 2), the initial value problem (2.1)-(2.2) has a unique local-in-time solution 
4> € (7 0 ([0, T); H 1 ( 


(2.3) 


piHonlf6lfTnTl[37|. 

NLS / GP is a Hamiltonian system, expressible in the form 


idt$> = 




with Hamiltonian energy functional 


H[F,F*]= / VF(x) ■ X7F*(x) + V(x) F(x)F*(x) + 
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(. F(x)F*(x)) a+1 dx. 


cr + 1 

The solution of the initial value problem satisfies the conservation laws: 

= H[$o] 

AA[$(-,t)] = AA[d> 0 ] 


where %[•] is defined in (2.4) and 

A f[F,F*] = [ F{x)F*(x)dx. 


(2.4) 


(2.5) 

( 2 . 6 ) 


(2.7) 


The conserved integrals Ji and M are associated, respectively, with the invariances t t+to, to £ M 
(time-translation) and 4> i-)- &e ie , 0£l. 

For the nonlinear Schrodinger equation (NLS), where V = 0, we have the additional symmetries 

1. 4 >(x,t) 4>(x + xq, t) , xo E (Translation invariance) 

2. 4>(x,t) i ^ [4>] (x,t) = 2£t,t), £ £ M d (Galilean invariance) 


For NLS with the homogeneous power nonlinearity in (2.1) the equation is also dilation invariant: 

<3?(x, t) i —> A° r< I>(Ax, A 2 t) (2-8) 
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The case where g < 0 is called the focusing case and the case where g > 0 is called the defocusing 
case. This terminology refers to the nonlinear term’s tendency to focus (localize) energy or defocus 
(spread) energy. Independently of the nonlinear term, the linear potential may be an attractive 
potential, which concentrates energy, or a repulsive potential which allows energy to spread. 

Concerning global-in-time behavior of solutions, the initial value problem is globally well-posed 
for any H 1 initial condition if either g > 0 or if g < 0 and a < 2/d (subcritical case). For g < 0 
and a > 2/d solutions with H l initial data may develop singularities in finite time (blow up). See, 
for example, 

3 Bound States - Linear and Nonlinear 

3.1 Linear bound states 

We start with the linear Schrodinger equation: 

idt* = (-A + V{x)) V = H (3.1) 

Here V denotes a smooth and real-valued potential satisfying 


V{x) —>• 0 sufficiently rapidly as |x| -» oo . (3-2) 

Bound states are solutions of the form T = e~ lEt if(x), where if satisfies the eigenvalue problem 

Hif = {-A + V)if = Eif, ifeH 1 ^). (3.3) 


We expect the eigenvalue problem ( |3.3| ) to have a non-trivial solutions if V is an attractive 
“potential well”; intuitively, the set where V(x) < 0 should be sufficiently wide and deep. 

The following result gives a condition on the existence of a ground state, a nontrivial solution 


of (3.3) of minimial energy, E H3CS1 . 


Theorem 3.1 Denote by 

Eo* = inf (Hf,f)= inf [ |V/(x)| 2 + V(x)\f(x)\ 2 dx. 

f R d l/l 2 =l f R d I/I =1 JM. d 

If —oo < < 0 then the infimum is attained at a positive function V’o* £ H 1 (M rf ) ; which is a 

solution of 

Hifo* = E 0ir if 0ir , E 0 * = (Hif 0 *,ifo*) 


3.2 Nonlinear bound states 

Here we consider spatially localized solutions of NLS/GP of the form 

T(x,t) = e~ lEt if E {x) , (3.4) 
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where ip E satisfies the nonlinear elliptic equation 


(-A + V)fi E + g\^E\ 2a ^E = Ei/j e , 


(3.5) 


Solutions of the nonlinear eigenvalue problem (3.5) are often called nonlinear bound states or 


solitary waves. A positive and decaying solution is often called a nonlinear ground state. In the 
case where V(x) is non-trivial, such states are often called nonlinear defect modes. In this case, 
the spatially varying potential, V(x), is viewed as a localized defect in a homogeneous (translation- 
invariant) background. 

The question of existence of nonlinear bound states of NLS / GP has been studied extensively by 
variational and bifurcation methods. See, for example, [991 151 IT081 1771 188] as well as the discussion 


in sections 3.6 m and 1 3.7 1 m- 

The following result, which we shall refer to later in this article, concerns bifurcation of a 
nonlinear ground state from the ground eigenstate of the linear operator — A +V; see, for example, 

mm 


Theorem 3.2 Assume V is real-valued, smooth and rapidly decreasing at infinity and that —A + V 
has at least one discrete eigenvalue. Let Eq*) denote the ground state eigenpair given by 


Theorem 3.T here, ipo* may be chosen to be positive and normalized in L 2 . Consider the equation 
for nonlinear bound states of NLS / GP with g = — 1 (attractive nonlinear potential): 

( -A + V(x) - \if E \ 2a ) fi>E = Efi E , E< 0, a > 1. (3.6) 

Then, there is a constant <5q > 0 and a non-empty interval T = [Eq* — ho,Eo,*) such that for any 


E G I, equation (3.6) has a positive solution, which for E tending to Eg* has the expansion: 


iPe(x) = p{E) ( ipo*(x) + 0(p(E) 2 *)), p(E) = |E 0 * —E|5 


3.3 Orbital stability of nonlinear bound states 


4 * 


2 ( 7 + 2 ) 


(3.7) 


In this section we discuss the orbital Lyapunov stability of nonlinear grounds states. In particular, 
we prove nonlinear stability in E 1 (M n ) modulo the natural group of symmetries. We shall in this 
section consider the case where g = — 1, the case of focusing nonlinearity. Given a positive H 1 


solution of (3.6) introduce its orbit: 


Ogs = {^ E (x)e^ : 7 g [ 0 , 2 tt )}, V^O 

O gs = {^ E (x + x 0 )e^ : T e [0, 2vr), x 0 G M n } , V = 0 . 


(3.8) 


We say that the ground state is orbitally stable’ in E 1 (M n ) if the following holds. If <I>(x, t = 0) 
is close to element of O gs in H l then for all t / 0, ^>(a:, t) is H 1 close to some (typically t dependent) 
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element of O gs . In order to make this precise, we introduce a metric which measures the distance 
from an arbitrary H 1 function to the ground state orbit: 


dist (u, O gs ) = inf ||w — 'I'oe* 7 ||#i, V ^ 0 
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dist (u, Og S ) = inf ||it(- + xq) — d'oe* 7 ||#i, V = 0. (3-9) 

7> x 0 

A more precise statement of H 1 — Lyapunov stability is as follows. For any e > 0 there is a <5 > 0 
such that if 

dist ($(-,0), Og S ) < 6 (3.10) 


then for all t ^ 0 

dist (<F(-,f), e> gs ) < e. 

Nonlinear bound states are critical points of the Hamiltonian T~i subject to fixed L 2 norm, AT. 
In particular, a nonlinear bound state, iJje, of frequency E satisfies S£e[/, = 0, where 

£ E [f} = H[f\ - E j |/| 2 

W\ = J |V/| 2 + H(x)|/| 2 -^ t |/| 2ct+2 

For subcritical nonlinearities, a < 2/d, stable nonlinear ground states may be realized as constrained 
global minimizers of the variational problem: 


/] subject to fixed A f[f] > 0. (3-11) 

Orbital H 1 stability is a consequence of the compactness properties of arbitrary minimizing se¬ 
quences mm 

However in general, depending on the details of the nonlinear terms, stable solitary waves may 
arise as local minimizers of % subject to fixed AT. We now discuss stability, in this more general 
setting. Introduce the linearized operators L±, which are real and imaginary parts of the second 
variation of the energy £e- 


L+ = - A + V(x) - (2(7 + 1)^ - E. 
L_ = -A + V(x) - - E 


(3.12) 

(3.13) 


We define the index of ^e by 


index (i/je) = number of strictly negative eigenvalues of L + . (3-14) 

Theorem 3.3 ( H 1 Orbital Stability) 1. Assume the conditions 
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(51) Spectral condition: i(e > 0 and index ('i)e) = 1 

(52) Slope condition 

-^U[^ E ] = J \if E (x)\ 2 dx < 0 (3.15) 

Then, ip E is H l orbitally stable. 


2. Assume i/j e > 0. If index (Vie) >2 or > 0 then i(e is linearly exponentially 

unstable. That is, the linearized evolution equation (see (3.19) below) has a spatially localized 
solution which grows exponentially with t. 


See [TUDj TIP , I5T11521 5D, T[1 • In |10t>| it was shown that d E A /"[V'e] > 0 implies the existence of an 
exponentially growing mode of the linearized evolution equation. 


We give the idea of the proof of stability. For simplicity, suppose V ( x ) is non-trivial, In this 


case the ground state orbit consists of all phase translates of i(e', see the first line of (3.8). Let e 


be an arbitrary positive number. We have for t ^ 0, by choosing 6 in (3.10) sufficiently small 

e 2 ~ £e[\I/(-,0)] - £e[iPe] 

= £ E [&(-,t)] — £e[^e], by conservation laws, 

= £ E [^(-,t)e n ] — £e['i(e\, by phase invariance 

= £e[i(e + u(-,t) + iv(-,t)\ — £e[^e] (definition of the perturbation u + iv, «,»eK) 

~ (L + u(t),u(t)} + {L-v(t),v(t)) (by Taylor expansion and S£e[^e] = 0) (3.16) 


If L + and L_ were positive definite operators, implying the existence of positive constants C+ and 
C- such that 


(L + u,u) > C + ||u||^i, 
(L-v,v) > C- \\v\\ 2 h i 


(3.17) 

(3.18) 


for all u, v £ H ] , then it would follow from (3.16) that the perturbation about the ground state, 
u(x,t ) + iv(x,t), would remain of order e in H l for all time 1 / 0. The situation is however 
considerably more complicated. The relevant facts to note are as follows. 


(1) L-ipE = 0, with i)e > 0. Hence, V’E is the ground state of L_, 0 G a(L-), and L_ is a 
non-negative with continuous spectrum [ |i7|,oo). 

(2) For small L 2 nonlinear ground states, L + has exactly one strictly negative eigenvalue and 
continuous spectrum [ |I7|,oo). 

The zero eigenvalue of L_ and the negative eigenvalue of L + constitute two bad directions, which 
are treated as follows, noting that u(-,t) and v(-,t) are not arbitrary H 1 functions but are rather 
constrained by the dynamics of NLS. 
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To control L_, we choose 7 (t) so as to minimize the distance of the solution to the ground state 


orbit, (3.9). This yields the codimension one contraint on v: (v(-, t), i\)e) = 0, subject to which 


(3.18) holds with C- > 0. 

To control L + , we observe that since the L 2 norm is invariant for solutions, we have the codi¬ 
mension constraint on u : (u(-, t), 4 >e) = 0. Although VtE is not the ground state of L +1 it can be 


shown by constrained variational analysis that if the slope condition (3.15) holds, then constraint 


on u places u in the positive cone of L + , i.e. C + > 0 in (3.17). Thus, positivity (coercivity) esti¬ 


mates (3.17) and (3.18) hold and £e serves as a Lyapunov functional which controls the distance 
of the solution to the ground state orbit. The detailed argument is presented in [ llO j. 


We also note that the role played by the second variation of £ in the linearized time-dynamics 
m\- Let \I >(t ) = {ipE + u + iv)e lEt . The linearized perturbation, (u(t), v(t)Y, satisfies the linear 
Hamiltonian system: 



with conserved (time-invariant) energy 


(3.19) 


Q(u,v) = (L + u,u) + (L-v,v) 


(3.20) 


The above constrained variational analysis underlying nonlinear orbital theory corresponds to the 
H 1 boundedness of the linearized flow (3.19), restricted on a finite codimension subspace. This 
subspace is expressible in terms of symplectic orthogonality to first order generators of symmetries 
acting on the ground state. This structure is used centrally in many works on nonlinear asymptotic 
stability and nonlinear scattering theory of solitary waves; see section [ 6 ] and references cited. 


In the following subsections, we give several examples: a) the free NLS soliton, b) the nonlinear 
defect mode of a simple potential well, c) nonlinear bound states and symmetry breaking for the 
double well, and d) gap solitons of NLS / GP with a periodic potential. 


3.4 The free soliton of focusing NLS: V = 0 and g — — 1 

In this case, there is a unique (up to spatial translation) positive and symmetric solution, a ground 
state, which is smooth and decays exponentially as |x| —> 00 . For results on existence, symmetry 
and uniqueness of the ground state see, for example, mumm- 

For example, the focusing one-dimensional cubic (a = 1) nonlinear Schrodinger equation 

id t ip = -d 2 x - IV’I 2 ^, (3.21) 


has the solitary standing wave: 


' i/j so i(x,t ) = e jt \/2 sech (x) , 


(3.22) 













By the symmetries of NLS (section [ 2 ]) , we have the extended family of solitary traveling waves 
of arbitrary negative frequency E = — A 2 , A > 0: 

G\,x 0 ,v,eH>soi\(x,t) = \e lX2t V2 sech (A (x - x 0 - 2vt)) e lv ( x ~ x o-vt) e i6.^ ^ ^ ^ e jj. (3.23) 


Theorem 3.3 implies that the positive solitary standing wave of translation invariant NLS is 
orbitally stable if a < 2/d and is unstable if er > 2/d PZH H3H ED- The solid curve in Figure 
3.5 shows the family of solitary waves, (3.23), of (3.21[) bifurcating from the zero solution at the 


edge-energy of the continuous spectrum of —d 2 . 


3.5 V(x), a simple potential well; model of a pinned nonlinear defect mode 

hi [55] the bifurcation of nonlinear bound states of NLS/GP from localized eigenstates of the linear 
Schrodinger operator — A + V was studied. The simplest case of such a nonlinear defect mode is 
for NLS / GP with V ( x ) is taken to be a Dirac delta function potential well: 


id t if = H'if — \i/)\ 2 if), H = — d'l - j(x), 7 > 0 . 


(3.24) 


For the well-posedness theory of the intial-value-problem for (3.24) in C'°([0, T): iL 1 (M)) see [60 


and, more specifically, section 8 D of m- In this case, the nonlinear defect mode of frequency 
E = —A 2 < —( 7 / 2) 2 is explicitly given by: 


r if so i(x,t-, 7 ) = \e lX2t V2 sech ^A|x| + tanh 1 ^0 


jo 


This family of nonlinear bound states, which is pinned to the “defect” at x = 0, bifurcates from 
the zero solution at E = E+ = — ( 7 / 2 ) 2 , the unique negative eigenvalue of iL; see the dashed curve 


in Figure 3.5 



E 


Figure 1: Squared L 2 norm, Af, vs. frequency for the free soliton (solid) and pinned soliton (dashed). 


9 















3.6 NLS / GP: Double well potential with separation, L 

Consider now the NLS / GP with a double-well potential, Vl(x): 

idt= —Ai/j — Vl(x)iJj — \ip\ 2, ip . 


We take the double-well, Vl(x) to be constructed by centering two identical single-bound state 


potentials, Vq(x ), of the sort considered in the example of section 3.5, at the positions x = ±L: 


V L (x) = V 0 (x + L) + V 0 (x-L) 


(3.25) 


A sketch of a one-dimensional double well potential, Vl(x), and the associated spectrum of H = 
—d 2 + Vl(x ) is displayed in the right panels of Figure [ 2 ] 



Figure 2: Left: Single-well potential and its spectrum: one discrete eigenvalue, marked by “o”, and continuous 
spectrum, R+. Right: Double-well potential and its spectrum: two nearby discrete eigenvalues (marked and “x”) 
and continuous spectrum, R+ 

As in the example of section [375] (see [55]) there are branches of nonlinear bound states which 
bifurcate from the zero solution at the discrete eigenvalue energies. We focus on the branch ema¬ 
nating from the zero solution at the ground state eigenvalue of Hl. For small squared L 2 norm, 
A f, this is the unique (up to the symmetry if; 1 —y e ld ip) nontrivial solution branch. This solution 
has the same symmetries as the ground state of the linear double well potential [52]. That is, for 
M small, ipE(jV) is bi-modal with peaks centered at x ± L. Increasing the L 2 norm (or its square, 
J\f) we find that there is a critical value J\f cr > 0, such that for A f > J\T cr there are multiple non¬ 
linear bound states; see Figure [3j In particular, at {E(J\f CT ),ijjj^ ci ) there is a symmetry breaking 


10 













bifurcation. Specifically, for J\f > A/" C r> there are three families of solutions: the continuation of 
the symmetric branch (dashed curve continuation of the symmetric branch) and two branches of 
asymmetric states, corresponding to solutions whose mass is concentrated on the left or right side 
wells. The bifurcation diagram in Figure [3] show only two branches beyond the bifurcation point. 
The solid leftward branching curve represents both asymmetric branches, one set of states being 
obtained from the other via a reflection about x = 0. 

Figure [3] also encodes stability properties. For 0 < AT < J\f cr , the (symmetric) ground state is 
stable, while for A f > A f cr the symmetric state is unstable, at A f = A/" cr stability is transferred to 
the asymmetric branches (solid curve for A f > A/" cr ). 



Figure 3: Bifurcations from the two discrete eigenvalues of the double-well: Bifurcation curve emanating from the 
zero state at the ground state (lowest) energy shows symmetry breaking at some positive Af cr - Bifurcation curve 
emanating from the zero state at the excited state energy shows no symmetry breaking. 

The above discussion is summarized in the following |65l [72]: 


Theorem 3.4 (Symmetry Breaking Bifurcation) Consider the nonlinear eigenvalue problem 

{-A + V L )f; E + g\if E \ 2 ^ E = E^ e , ip E € H 1 (M d ) (3.26) 


where g < 0 and Vl denotes a double well potential with separation paramater, L, as in (3.25). 
Denote the ground and excited state eigenvalues of —A + Vl by: 


Fq*(L) < Eii'(L) < 0. 


(3.27) 


LetM[f] = f Rd |/| 2 . 

There exists a positive constant, Lq, such that the following holds. For all L > Lq, there exists 
A/cr = A/" cr (T) > 0 such that 

(a) For any M < A f cr there is a unique non-trivial symmetric state. 
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(b) The point (EMe) = (E CT ,ipE cr ) * s a bifurcation point i.e. 

there are, for AT > M cr , two bifurcating branches asymmetric states, consisting respectively of 
states concentrated about x = EL. 


(c) Concerning the stability of these branches: 

(cl) M < M cr : The symmetric branch is orbitally stable. 

(c2) M > M cr '- The symmetric branch is linearly exponentially unstable; 


these states have index = 2; see (3.14). 


(c3) M > M cr : The asymmetric branch is orbitally stable, 
(d) Symmetry breaking threshold: 

M cr (L) ~ E lit (L) - Eo*(L) 


(3.28) 


(The gap, Eu(L) — E$+{L) is exponentially small for large L 


'■) 



Figure 4: Symmetric, asymmetric and anti-symmetric modes arising in the bifurcation diagram of Figure [ 3 ] 


Theorem 3.4 is proved by methods of bifurcation theory 


Specifically, for L sufficiently 
large we can use a Lyapunov-Schmidt reduction strategy to reduce the nonlinear eigenvalue problem 
to a weakly perturbed system of two nonlinear algebraic equations depending on M. The unknowns 
of this system are essentially the projections of the sought nonlinear bound state onto the linear 
ground and excited states of Hi. The bifurcation structure of this pair of nonlinear algebraic 
equations is of the type displayed in Figure [3| This bifurcation structure is then shown to persist 
under the (infinite dimensional) perturbing terms to this finite algebraic reduction. The latter are 
controlled via PDE estimates and an implicit function theorem argument. 

An illustrative study of the global bifurcation structure for the exactly solvable case where 
Vl{x) is taken to be a sum of two attractive Dirac-delta wells centered at x = ±L is presented 
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in [60]. In [[63] the local bifurcation methods described above are used to study the NLS / GP 
with a triple well potential. A recent study of global properties of the bifurcating branches and 
their stability properties for general nonlinearities, |'0| 2cr '0, with a > 1/2 is presented in [65J. In 
this work the authors prove, for the case of symmetric double-well potentials, that as L varies a 
symmetry breaking bifurcation occurs once the potential develops a local maximum at the origin. 
The analysis in [65] allows for the bifurcation to occur at large J\f cr , outside the weakly nonlinear 
regime where the analysis in m applies. 

Local bifurcation methods can be applied as well to study excited state branches. An example 
is the branch emanating from the linear excited state energy, along which it can be shown, for small 
AT, that there are no secondary bifurcations; see Figure [3] and the corresponding mode curve which 
changes sign in Figure [4] 

An approach which complements the local bifurcation approach is variational. For example, in 
[2| symmetry breaking in a three-dimensional nonlinear Hartree model having an attractive non¬ 
local potential is proved for A f sufficiently large. Their arguments are straightforwardly adaptable 


to the variational formulation (3.11) and yield that for a < 2/d the ground state is asymmetric. 


The proof is based on a trial function argument and the intuition that for A f sufficiently large it is 
energetically preferable to concentrate mass in the left well or in the right well but not both. 


3.7 NLS / GP: V(x) periodic and the bifurcations from the spectral band edge 

We now consider the nonlinear bounds states of NLS / GP for the case where V is periodic on M rf . 
In this case, the spectrum of H = —A + V(x) is absolutely continuous. Let Vv denote a nonlinear 


bound state of frequency fj,. Any bifurcation of nontrivial solutions of (3.5) from the zero solution 


must occur at an energy in the continuous spectrum. Indeed, the example of the free soliton of 
section [3A| is an example; V = 0 is periodic (!) and has continuous spectrum equal to the half-line 
M_l_ = [0, oo). The soliton bifurcates from the zero-amplitude solution (as measured say in L°° or 
L 2 ) from the continuous spectral edge (solid curve in Figure [ 3 ]). Indeed, for NLS with V = 0 and 
power nonlinearity — \'iJj\ 2rT ip 

= (-yu)^Y’-i(\/ = 7^) (3.29) 


and therefore 


A = 


Mi* 


= (-/lK 




(3.30) 


For V(x) a general periodic potential the spectrum of H = — A+V is the union of spectral bands 
obtained as follows. Consider the d— parameter family of periodic elliptic eigenvalue problems: 


^ — (V-Hfc) 2 + V(x)^ju(x) = n u(x) 
u(x) is periodic with the periodicity of V(x) , 


(3.31) 
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where k varies over a fundamental dual lattice cell, the Brillouin zone, B C M rf . For each fixed 
k G B , the spectrum of (3.31|) is discrete and consists of eigenvalues denoted: 


Pi(k) < < ■ ■ ■ < M&(fc) < ■ ■ • , 


(3.32) 


listed with multiplicities and tending to positive infinity. The corresponding eigenfunctions are 
denoted, pj(x ; &). The states {e lkx pj(x; k)} where j > 1 and k & B are complete in L 2 (M rf ). 

It can be shown that nonlinear bound states bifurcate from edges of the spectral bands: to the 
left, if the nonlinearity is attractive (g < 0 ) and to the right if the nonlinearity is repulsive (g > 0 ) 
153 EH ESI- We expect similar scaling of the L 2 norm of such more general edge-bifurcations: 

~ \P>~P*\"~* (3-33) 

where \p — p+\ is the distance to the spectral band edge located at p*. Now the details of the 
bifurcation at the edge depend on the periodic structure. We first give a heuristic picture and then 
state a precise theorem. 

For \p — p+\ small, that is for p near a spectral band edge, the nonlinear bound state, x ), 
should be exponentially localized with decay ipn(x) ~ exp (—\p — p+\z\x\) for |x| —> oo. Due to 
the separation of length scales: |p — p+ 1 _1 S> period of V, we expect 'ip^(x) should oscillate like 
the Floquet-Bloch mode pi(x;0) = p*(x), associated with the band edge energy, //*, with a slowly 
varying and spatially localized amplitude, F(x), : 

5*F(5x) x p*(x), 5 = \p — p*\^ , (3.34) 


where F(x) satisfies an effective medium (homogenized / constant coefficient) NLS equation. 

To simplify the discussion we will focus on bifurcation from the lowest energy (bottom) of the 
continuous spectrum, p+ = pi(0), for the case of an attractive nonlinearity (g < 0). Before stating 
a precise result, we introduce the inverse effective mass matrix associated with the bottom of the 
continuous spectrum defined by 

d 2 pi(k = 0)\ 

dkidkj ) ^ ^ 

4 {d Xj p* , (—A + V - p+y 1 d Xi p±) 

(P* , P*) 


(3.35) 


A eff = ^D 2 p 1 (k = 0) = ^ 


and F e ff, the effective nonlinearity coefficient defined by 

/ p*(x ) 2cr+2 dx 


F e ff — 


J p*(x ) 2 dx 


> 0 . 


(3.36) 


The matrix A e g, is clearly symmetric and, for the lowest band edge, it is positive definite m- 
The following result, proved in 153 , describes the bifurcation from the bottom of continuous 
spectrum (left end point of the first spectral band). See also [54, [§B1,89]. A proof of the case of a 
general band edge is discussed in [59]. 
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Theorem 3.5 (Edge bifurcations of nonlinear bound states for periodic potentials) Let 

V{x) denote a smooth and even periodic potential in dimension d = 1,2 or d = 3. Let x$ denote 
any local minimum or maximum ofV(x). Denote by /x* = inf<r(—A +V). Let Fa,y{u) the unique, 
centered at y = 0, positive solution (ground state) of the effective / homogenized nonlinear 

Schrodinger equation with inverse effective mass matrix, A, and effective nonlinearity, F. 

- £ ir i A ” 1 «ir F(y) - r *» ■ < 3 - 37 ) 

ij =1 1 3 

Then, there exists 5 q > 0 such that for all /x 6 (/x* — <fo, /x*) there is a family of nonlinear bound 
states p t-A V’m( x ) 

- (p*-p)^ F AeS} r eft ( VM* - T (x - x 0 ) ) p*(x) ->• 0, as /x f in H 1 (R d ) . (3.38) 


The proof of Theorem 3.5 [59], like the study of bifurcations from discrete eigenvalues, proceeds 
via Lyapunov-Schmidt reduction and application of the implicit function theorem. However, while 
the analysis of bifurcation from discrete spectrum leads to a finite dimensional (nonlinear algebraic) 
bifurcation equation, bifurcation from the continuous spectrum leads to an infinite dimensional 


bifurcation equation, a nonlinear homogenized partial differential equation (3.37); see also 
Examples of recent applications of these and related ideas to other systems appears in [291 ESI EH 
[MU35U3UE2U33]. 


Since (3.37) is a constant coefficient equation, F A ,v can be related, via scaling, to the nonlinear 
ground state of the translation invariant nonlinear Schrodinger equation. Thus, the shape of the 
bifurcation diagram /x i-)- A f[ipn\ (see Figure [3]) can be deduced for /x near /x*; see [59] . We also note 
that the stability / instability properties of ground states of (3.37) for the time-dependent effective 
nonlinear Schr’odinger equation, are a consequence of Theorem , 


3.3 


in the translation invariant case. 


4 Soliton / Defect Interactions 

In this section we turn to the detailed time dynamics of a soliton-like nonlinear bound state in¬ 
teracting with a potential. This question has fundamental and applications interest. From the 
fundamental perspective it is an important problem in the direction of developing a nonlinear 
scattering theory for non-integrable Hamiltonian PDEs. And, from an applications perspective, 
nonlinear waves often arise in systems with impurities, e.g. random defects in fabrication or those 
deliberately inserted into the medium to influence the propagation; see, for example, [35], 

Referring to Figure [4] we sketch the different stages present in the dynamics of a soliton which 
is initially incident upon a potential well: 

The results of formal asymptotic and numerical studies of soliton / defect interactions (see, for 
example, [MUMUn]) suggests a description of the time-dynamics in terms of overlapping time 
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Figure 5: P = A f[i> fj] vs. g for the Id - NLS / GP with periodic potential, V(x ) for the case a = cr cr = 2/d = 2. 
Dashed line P = P cr is the (g— independent) value of Affi/v] for the case V = 0. Solid curve is the curve g H > A/"[VvL 
where g < g* = inf(cr(—A + F)) for the case where F is a non-trivial periodic potential. Nontrivial effective mass 
matrix implies P+ = lim^* A/"[VV»] < P C r- See [59] . 


epochs. Very roughly, these are: 

(i) short and long / intermediate time-scale classical particle-like dynamics, which govern motion 
and deformations of the soliton as it interacts with the defect and exchanges energy with the defect’s 
internal modes (pinned nonlinear defect modes) and 

(ii) the very long time scale, during which the system’s energy resolves into outgoing solitons moving 
away from the defect, small amplitude waves which disperse to infinity and an asymptotically stable 
nonlinear ground state of supported in the defect. A key process in this asymptotic resolution is 
radiation damping due to coupling of discrete and continuum degrees of freedom. See sections [5] 
and El 

Examples of rigorous analyses for transient / intermediate time-scale regimes: 


(i) Soliton scattering from a potential well: Detailed reflection, transmission and trapping for 
solitons incident on a potential barrier or potential well mmmm- 


(ii) Soliton evolving in a potential well: In m the evolution of an order one soliton in a potential 
well is considered. The detailed nonlinear breathing dynamics, as the soliton relaxes toward 
its asymptotic state are considered in [58]; see also section [6] and the related work [351 HOI ITT] , 


In (80j the evolution of a weakly nonlinear (small amplitude) NLS / GP solutions in a double¬ 
well, for which there is symmetry breaking (see section 3.6) is studied. Figure [4] displays 
phase portrait of the reduced Hamiltonian dynamical system in [801 . Orbits around the left 
(respectively, right) equilibrium map to a soliton executing a long-time nearly periodic (back 
and forth) motion with the left (respectively, right) well of the double well. This analysis has 
been recently extended to certain orbits outside the separatrix m- 
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Position 


Figure 6: A soliton interacts with a defect, support around x = 0. As time advances, the coherent structure 
interacts with the defect. Part of the incident energy is scattered to infinity and part of it is trapped in the defect. 
The trapped energy settles down to a stable nonlinear defect mode. 

For 1 > 1 the center of mass eventually crystallizes on a stable nonlinear ground state. For 
A f > N cr this asymptotic state will be an asymmetric state (section [376] ) centered in the left or 
in the right well. See Figure [4] which shows the computed center of mass motion for a solution 
of NLS / GP; initially there is oscillatory motion among the left and right wells. However, 
during each cycle some of the soliton’s energy radiated away to infinity. The corresponding 
motion in the reduced phase portrait Figure [4] is damped and the actual center of mass 
trajectory is transverse to the level energy curves. Eventually the separatrix is crossed and 
the solution settles down to a stable asymmetric state | T9l SUET] We emphasize that this 
picture is only heuristic although there has been considerable progress toward understanding 
the radiation damping and ground state selection in such systems. We now turn to these 
phenomena in the next section. We conclude this section noting work on the dynamics of 
forced / damped NLS and its reduced finite dimensional dynamics, capturing regular and 
chaotic behavior; see, for example, mm- 

5 Resonance, radiation damping and infinite time dynamics 

Our next goal is to discuss results concerning the infinite time dynamics of NLS / GP. We begin 
with an informal discussion of the linear Schrodinger equation (NLS / GP with g = 0): 

idti> = (—A + V(x))^ (5.1) 

Assume V(x) is a smooth and sufficiently decaying real-valued potentia for which —A + V has 
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Figure 7: Periodic dynamics of center of mass dynamics in the reduced (approximate) finite-dimenstional Hamil¬ 
tonian system [80]. Equilibria corresponding to stable asymmetric states centered on left and right local minima of 
double well. Oscillatory orbits decay to stable equilibria for the full, infinite-dimenstional NLS/GP. 


Plot of the Center of Mass for a solution u versus time t. 



Figure 8: Numerically computed dynamics of center of mass in the NLS evolution [SD]. Corresponding trajectory 
in the reduced phase portrait of Figure [f] cycles around several times outside the separatrix, crosses the separatrix 
and then slowly spirals in toward the right stable equilibrium. 
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bound states : 


(-A + V)i(>j = Sljifij, j = 0,1 ,..., il>j € L 2 
ip(x,t) = e~ lCljt tpj(x) 


Then, for sufficiently localized initial conditions, -0(0), the solution of the initial value problem 


for (5.1) may be written: 


= e * ( A+V ) Vo = %Qjt 'ip j (x)+ R(x, t) 


(5.2) 


The sum in (5.2) is quasi-periodic while the second term, R(x,t ), radiates to zero in the sense 
that for some a > 0 

II LP(R d ) < as i-> Too (5.3) 


Thus, as t —> oo, the solution ip(x,t ) tends to an asymptotic state which is quasi-periodic in time 
and localized in space. 


Question: What is the large time behavior of solutions of the corresponding initial value problem 
for the nonlinear Schrodinger / Gross-Pitaevskii, g ^ 0? 

For small amplitude it is natural to attempt expansion of solutions in the basis of localized 
eigenstates and continuum modes, associated with the unperturbed (solvable) linear Schrodinger 
equation. In terms of these coordinates, the original Hamiltonian PDE may be written as an 
equivalent dynamical system comprised of two weakly coupled subsystems: 


• a finite or infinite dimensional subsystem with discrete degrees of freedom ( “oscillators” ) 

• an infinite dimensional system (wave equation) governing a continuum “field” . 


These two systems are coupled due to weak nonlinearity. If one “turns off” the nonlinear 
coupling there are localized in space and time-periodic solutions, corresponding the eigenstates of 
the linear Schrodinger equation. If nonlinearity is present, new frequencies are generated, and these 
may lead to resonances among discrete modes or between discrete and continuum radiation modes. 
The latter type of resonance plays a key role in understanding the radiation damping mechanism 
central to asymptotic relaxation of solutions to NLS / GP as t ±oo and, in particular the 
phenomena of ground state selection and energy equipartition. 

These results will be discussed in section [6l A variant was studied in the context of the nonlinear 
Klein Gordon equation in |951 4]. In the following two subsections we present a very simple example 
of emergent effective damping in an infinite dimensional Hamiltonian problem. 
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5.1 Simple model - part 1: Resonant energy exchange between an oscillator 
and a wave-field 

We consider a solvable toy model of an infinite dimensional Hamiltonian system comprised of two 
subsystems, one governing discrete and the other governing continuum degrees of freedom. We shall 
see that coupling of these subsystems leads to energy transfer from the discrete to the continuum 
modes and the emergence of effective damping. Our model is a variation on the work of Lamb 
(1900) [75] and Weisskopf-Wigner (1930) [ 111] : see also [93] HH ESI mi EZl [Ml EM HI] 

Consider a system which couples an “oscillator” with amplitude a(t) to a “field” with amplitude 
u(x, t), f G R, x € M: 

da £ {t) 


dt 


+ iua £ (t) = — eu £ (0, t), 


dtu £ (x,t) + cd x u £ (x,t ) = e5(x)a £ (t). 


(5.4) 

(5.5) 


Here, e and c are taken to be a real parameters, say with c > 0, and 5(x) denotes the Dirac delta 
function. The amplitudes a(t) and u(x, t ) are complex-valued. 

Note that the dynamical system (5.4)-( ph5] ) conserves an energy, dt£[a(t),u(-,t)] = 0. 

(5.6) 


£[a(t),u(-,t)} = |a(f)| 2 + / \u(x, t)\ 2 dx 

Jr 


We consider the initial value problem for (5.4)-(5.5) with initial data 

a(0) arbitrary and u(x, 0) = 0 ; 


(5.7) 


we perturb the oscillator initially, but not the field. 

For e = 0 the oscillator and wave-field are decoupled and there is an exact time-periodic and 
finite energy global solution, a bound state: 


,£= 0 ^ 

,£= 0 / 


r ”(x,f) 


= e 


—iujt 


o 

a(0) 


a £=0 (f) 

For £ / 0, the oscillator and field are coupled; as time evolves energy can be transfer among 
the discrete and continuum degrees of freedom. It is simple to solve the initial value problem 
for any choice of e and arbitrary initial data, e.g. by Laplace transform. In particular, one can 


proceed by solving the wave equation, (5.5), for u £ = u(x,t]a £ ) as a functional of discrete degree 
of freedom a £ (t ) and then substitute the resulting expression for u £ (0,f; a £ (t)) into the oscillator 
equation (|5.4|). The result is the closed oscillator 


+ iua £ (t) = ~—a £ (t ), t > 0. 
dt 2c 


with an effective damping term, whose solution is: 


a £ {t ) = e~ iut e-^ a 0 = t > 0. 


(5.8) 


(5.9) 
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That is, in terms of the oscillator, the closed infinite dimensional conserved system for a £ (t) and 
u £ (x,t ) may equivalently be viewed in terms of a reduced open and damped finite dimensional 
system for a e (t). 

Solving ( |5.8| ) for a £ (t) and substituting into the expression for u £ = u(x,t',a £ ), one obtains: 

0, x > ct> 0 


u e (x,t ) = 


..-Mt-f) e -fc(i-t) x<ct 


To summarize, for e = 0, the decoupled system has a time-periodic finite energy bound state. 
For the coupled system, £ ^ 0, the bound state has a finite lifetime; it decays on a timescale of 
order e -2 . The bound state loses its energy to the continuum degrees of freedom; the lost energy 
is propagated to spatial infinity. 


5.2 Simple model - part 2: Resonance, Effective damping and Perturbations 
of Eigenvalues in Continuous Spectra 


It’s physically intuitive that an oscillator coupled to wave propagation in an infinite medium will 
damp. We now connect this damping to the classical notion of resonance. 

We write the oscillator / wave-field model in the form 

idtip = A(X) V’, (5.10) 


where 


Let 


V>(t) = 


u(x, t ) 
a(t) 


, A(e) = 


—icd x +ie5(x ) 
to 


m = e- i£ Vo = 

This yields the spectral problem 

A(e)^o = Eij ;o , 

or equivalently 

E 


U ( X A) \ —iEt ( u 0 
a(t) 1 


ao 


—icd x +ie5(x) 
This spectral problem may be considered on the Hilbert space: 


u 0 

ao 


u o 
« o 


H = 


u 


(E L 2 (M) x C : || (u, a)< oo > , 


where 


\\(u,a)\\ n = / \U(x)\ 2 dx + |a| 2 < oo 

Jr 

We note that for e = 0 (decoupling of oscillator and wave-field), we have 

E 

A{e = 0) is diagonal and has spectrum given by: 


Uq 

ao 


—icd x 0 
0 to 


u 0 

ao 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 
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a point eigenvalue at E = ui with corresponding eigenstate ij^ = 


continuous spectrum M = {E = ck : k E M} with eigenstates ipk = 


and 


0 ikx 


o 


k e 


Thus, A(e = 0) has an embedded eigenvalue, u>, in the continuous spectrum, M, and the source the 
damping is seen to be this resonance. Understanding the coupling of oscillator and wave-held is 
therefore related to the perturbation theory of an embedded (non-isolated) point in the spectrum; 
see, for example, [87, '53| [93] and references therein. 

In this present simple setting, this perturbation problem can be treated as follows. Note that 
for £ / 0, since c > 0 we expect that energy is emitted by the oscillator and it gets carried to 
positive infinity through u(x,t). Thus we expect that on a fixed compact set, u(x,t), will decay to 
zero. Alternatively, if we choose a > 0 and study U(x,t ) = e~ ax u(x,t), U(x,t) can be expected to 
decay to zero as t advances. Indeed, it can be checked that U(x,t) satisfies a dissipative PDE, for 
which the L 2 (M) norm of U(x,t) decays as t increases; see the discussion of section 5.1 


This behavior is reflected in the spectral problem. Consider the of change variables (u, a) i —> 


( e ax U,a ) in the spectral problem (5.14). This yields 


E 


U 0 

ao 


—icd x — ica +iee ax S(x ) 
—is (5{-)e a ', •} oj 


U 0 

ao 


= A a (e) 


U 0 

ao 


(5.18) 


which when considered on T 2 (M) is equivalent to the spectral problem (5.14) for ( u(x),a ) in the 
weighted function space: 


where 


For e = 0, 


U a = {(u(x),a) : \\u\\n a < oo },, 

|it||% a = / \u(x)\ 2 e~ 2ax dx + |a| 2 < oo , a > 0. 

Jr 


E 


vo 

ao 


—icdr — ica 0 


0 


u J \ a 0 


v 0 


Aq,(0 ) is diagonal and has spectrum: 

• point eigenvalue at E = u> with corresponding eigenstate = 


(5.19) 

(5.20) 

(5.21) 


continuous spectrum in the lower half plane along the horizontal line {E = ck — ica : k 6 

Akx 


with eigenstates V’fc = 


k e 


Note that for a > 0, u is an isolated eigenvalue of A a (e = 0) and we can therefore implement a 
standard perturbation theory to calculate the effect of e ^ 0 on this eigenvalue. 
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The first equation of (5.18) is 


(— icd x — ica) vq + ied(x)e ax ao = Ev o 


(5.22) 


This implies for i / 0 


v 0 (x) = 


e i( f x ~») x Vo+ x>0 

e i(f x -a)x Vo _ x < 0 


Integration of (5.18) over a small neighborhood of x = 0 yields 


The second equation of (5.18) implies 


^o+ — ^o- — +~ao 
c 


IE 


(E-u)ao = -—(v 0 ++v 0 -) 


Now choose uq+ / 0 and vq~ = 0. Then, 


(5.23) 


(5.24) 


E e = u - zT £ , T £ = — > 0 . 

2c 


The corresponding eigenstate is 


vb, a [x) = 


iui 7. P Hb- ac ) x : x > 0 
x < 0 


+fao e“ced2c 

0 , 


(5.25) 


(5.26) 


£ 

2c 1 ' 


which is in L 2 (M) provided a > 

For a = 0, we have 

«§(*) = 


+ ^a 0 e^cecS*, 


0 , 


x > 0 
x < 0 


The solution (5.27) satisfies an outgoing radiation condition at x = +oo and is in L 2 


(5.27) 
3 ~ 2ax dx) 


2 

provided a > . 

To summarize: the e = 0 eigenvalue problem has discrete eigenvalue, E° = u, embedded in the 
continuous spectrum, M with corresponding eigenstate f ^ j . This eigenvalue perturbs, for e / 0, 
to a complex energy in the lower half plane, E £ = u — iT e , T e > 0 with corresponding eigenstate 

( X ) \ 

J which solves the eigenvalue equation with outgoing radiation condition at x = +oo. 
The complex energy, E £ , may be viewed as a eigenvalue with normalizable eigenstate, vo(x) in the 
weighted function space, T~L a . Since E £ is in the lower half plane, the corresponding time-dependent 
state is exponentially decaying as t increases. 


Remark 5.1 Finally, we remark that the complex energy E £ may he viewed as a pole of the Green’s 
function, when analytically continued from the upper half E— plane to the lower half E— plane. 
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6 Ground state selection and energy equipartition in NLS / GP 


In this section we return to the question raised in section [5} here restated. For simplicity we consider 
NLS / GP for the case with cubic nonlinearity (<r = 1): 

id t <b = ( -A + F(x))<f> + 5 |$| 2 d>, i£l 3 (6.1) 

We assume that —A + V has multiple independent localized eigenstates. For simplicity we assume 
two distinct eigenvalues. 

As noted, the linear time-evolution (g = 0), for t 3> 1, settles down to a quasi-periodic state, 
a linear superposition of time-periodic and spatially localized solutions, corresponding to the in¬ 
dependent eigenstates. Note also, that for j / 0 there may be multiple co-existing branches of 
nonlinear defect states; see, for example, the discussion of bound states for the case where V is a 
double-well potential, discussed in section [3] 

Question: What is the long term (t f oo) behavior of the NLS / GP (g ^ 0), 
of small norm? 

Our results show that under reasonable conditions (which have been explored experimentally [78j, 
E3) we have: 

1. Ground state selection [!M, [57J, [73j/: The generic large time behavior of the initial value 

problem is periodic and is, in particular, a nonlinear ground state of the system. See also 
P 011011105119]. 

2. Energy equipartition \-j3tf : For initial data conditions whose nonlinear ground state and excited 
state components are equal in L 2 , the solution approaches a new nonlinear ground state, whose 
L 2 norm has gained an amount equal one-half that of the initial excited state energy. The 
other one-half of the excited state energy is radiated away. 

To state our results precisely requires some mathematical setup. For simplicity we assume that 
the linear operator —A + V has the following properties: 

(VI) V is real-valued and decays sufficiently rapidly, e.g. exponentially, as \x\ tends to infinity. 


(6.1) for initial data 


(V2) The linear operator —A + V has two eigenvalues £o* < Ed* < 0. -So* is the lowest eigenvalue 
with ground state V’o* > 0, the eigenvalue E i* is possibly degenerate with multiplicity N > 1 
and corresponding eigenvectors £i*,£ 2 *, ■ ■ -,Cn*- 

(V3) Resonant coupling assumption 

to* = 2-Si* — -Eq* > 0, . (6-2) 
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We remark on (V3). An important role is played by the mechanism resonant coupling of discrete 
and continuum modes and energy transfer from localized states to dispersive radiation. In sections 


5.1 and 5.2 we introduced and analyzed a simple model of this phenomenon. This model was linear 
and the resonance was due to a discrete eigenvalue, uj, embedded in the continuous spectrum. For 


NLS / GP, (6.1) resonant coupling of discrete and continuum modes arises due to higher time- 
harmonic generation by nonlinearity. The assumption (V3) is made to ensure that this coupling 
occurs at second order in (small) energy of the solution. However coupling at arbitrary higher order 
can be studied using the normal form ideas developed in (39J, 16, U ITT] . 


6.1 Linearization of NLS/GP about the ground state 


Recall Theorem 3.2 on the family of nonlinear bound states, E e -> ifE, bifurcating from the ground 


state. We now consider the linearized NLS / GP, (6.1), time-dynamics about this family of solutions. 
Let 

<h(x, t) = e~ lEt ( V’E +u+ iv ), 

where u and v are real and imaginary parts of the perturbation. Then the linearized perturbation 
equation about i/je, is 

!(:) 

where 


= JH{E) 


(6.3) 


L{E) = 


0 


L-{E) 


0 1 

-1 0 


L+(E) 0 
0 L_(E) 


= JH(E) 


(6.4) 


(6.5) 


-L+(E) 0 

The operators L + and L_ are given by: 

L-(E) = -A -E + V + g{^E? 

L+{E) = -A-E + V + 3g^ E ) 2 

The following result on the linearized matrix-operator, L(E), proved by standard perturbation 
theory m, is given in [42] (Propositions 4.1 and 5.1): 

Lemma 6.1 Assume (VI), (V2) and (V3). Let L(E) denote the linearized operator about a state 
on the branch of bifurcating bound states, given by Theorem \ 3. 

_ 1 

iPe(x) = p(E) ( i)o*(x) + O ( p(E) 2 ) ) , where p{E) = \ E 0 ± - E\* ^ \g\ J (6.6) 

For E = Eo+, the matrix operator L(Eq+) has complex conjugate eigenvalues ±z/3* = ±i(£i* — Eq+), 
each of multiplicity N. For \Eq* — E\ and small these perturb to (possibly degenerate) eigenvalues 
±i(3i(E ),..., ±i/3]\r(E) with corresponding neutral modes (eigenstates) 


6 

±*771 


6 

±im 


6 v 

±h]N 
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satisfying {f,m,Vn} = $m,n, (^m,fpE) = (Vm,d E 'lpE) = 0 . 

Moreover, 

0 f lim £ n = lim r] n G span{£i*. £ 2 *, ■ • •, 6v*} 

E — >Eo* E — >Eq+ 

in Sobolev H k spaces for any k > 0. 

Furthermore, we note that for \Eq* — E\ sufficiently small 


(6.7) 

( 6 . 8 ) 


2(3 n {E) + E ~ 2(l?i* — E$f) + Eq± — 2E\+ — Eq+ >0, n — 1, 2, • • •, N, 


(6.9) 


Remark 6.1 Equation (6.9) (see also (V3), (6.2)J ensures coupling of discrete to continuous spec¬ 


trum at second order in \E — E^\. 


6.2 Ground state selection and energy equipartition 

In this section we give a detailed description of the long term evolution. 


Theorem 6.1 (Ground State Selection) Consider NLS / GP, (6.1), with linear potential sat¬ 


isfying (VI), (V2) and (V3), and cubic nonlinearity (a = 1). Assume, most notably, that the non¬ 
negative (Fermi golden rule) expression for Tq{z,z*) in ( | 6. 15 ) is strictly positive. (See 1961 fj2l (fffj 
for detailed statement of all technical assumptions. This expression is always non-negative and 


generically strictly positive due to the assumption (6.2). ). 
Take initial conditions of the form: 


f>o(x) = e no [ if Eo + a 0 ■ £ + i/3 0 ■ g + R 0 ] 


( 6 . 10 ) 


where 70 and E$ 6l = (To* — So,Eo*) are real constants, «o and /3 q are real 1 x N vectors, and 
Ro : M 3 —> C, are such that 

\E 0 - E 0 f < 1, |Q "°| 2 , + ||^° |2 «1, \\(x) 4 Ro\\ H 2 < |a 0 | 2 + |A)| 2 (6.11) 

11 V’-Eo 112 

Then there exist smooth functions Eft ) : M + —> T, 7 ft) : M + —> M, zft) : M + —t C N and 
R{x, t) : M 3 x M + —» C such that the solution of NLS / GP evolves in the form: 

if{x,t) = g-i/o^V 7 ^) 

X [f>E{t) + Oi (z, z) d E ip E + ia 2 (z, z) if> E + {Re z[z, z]) • f + i ( Imz[z, z]) ■ g + R] , 

( 6 . 12 ) 

where lirn^oo Eft) = E^, for some E^ G T. 

Here, a±(z,z), a 2 {z,z ) : x C N —t M and z — z : C N x —> C N are polynomials of z and z, 

beginning with terms of order \z\ 2 . 
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(A) The dynamics of mass/energy transfer is captured by the following reduced dynamical system 
for the key modulating parameters, E(t) and z(t): 

\\^E(t)\\l = z*T 0 (z,z) z+ S E (t), (6.13) 

^ \z(t)\ 2 = -2z*T 0 (z,z) z + S z (t) . (6-14) 

Here, To is a Fermi golden rule (damping) matrix given 

T 0 {z,z*) = g 2 $ ([-A + V - uj+ - iO] -1 tf E (z ■ £) 2 , 4>e {z ■ £) 2 ) > c 2 \z\ 4 (6.15) 


and S E (t), S z (t ) < (1 + f) p , with / 0 °° S E (r)dT and / 0 °° S z {r)dT = o(|z 0 |) 2 . 

(B) Estimates on z(t ) and the correction R(t): For all t > 0, we have ||i?(t )||^2 < Coq. 
Moreover, the following decay estimates hold: 


(1 +x 2 )~ u R(t) 


< C(\\(x) 4 Mm) (l + t)~\ 


z(t )I < C (ll(x)Voll^) (1 + i)"5 


(6.16) 

(6.17) 


Theorem 6.2 ( Mass / Energy equipartition) Consider NLS / GP, (6.1), under the techni¬ 
cal hypotheses of Theorem 6.1 and assumptions on the initial data (6.11), where |ao| 2 + |/?o| 2 is a 
measure of the neutral modes ’ perturbation of the ground state. Recall that Eo is the energy of the 
initial nonlinear ground state and E^ is the asymptotic energy, guaranteed by Theorem \6.1\ Then, 
in the limit as t —>• oo, one half of the neutral modes’ mass contributes to forming a more massive 
asymptotic ground state and one half is radiated away as dispersive waves: 


1 


Too 11 2 — W^EoWl + z [ |ao| 2 + |A)| 2 ] + 0 ( |ao| 2 + |A)| 2 ) • 


(6.18) 


Figure [9] illustrates the phenomena of ground state selection and mass / energy equipartition in 
NLS / GP. These simulations are for the case of double-well potentials (see section 3.6). Asymptotic 
energy equipartition has been verified for the case when the stable ground state is symmetric 
(M < M C r) and the case when the stable ground states is asymmetric (A f > Af cr )- 

A result on equipartition of energy holds as well, for NLS-GP with a general power nonlinearity, 
(2.1), under more restrictive hypotheses [33]. Finally, we remark that it would be of interest to 


establish detailed results on energy-transfer in systems with multiple bound states for subcritical 


assumed to be strictly positive, is non-negative by the following: 

[—A + V — oj* — iOl” 1 = lim ( [—A + V — uo* — i<51 -1 — [—A + V — w* + ) 

2i «4.o v ’ 

= 7T S (—A + V — w*), where on G a con t (—A + V). 

Note: Sr [—A + V — w* — *0] _1 projects onto the generalized mode at energy w* € a CO nt (—A + V). 
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Figure 9: Numerical verification ground state selection and asymptotic energy equipartition. Computations 
performed by E. Shlizerman . 


nonlinearities. In this setting the current perturbative treatment of small amplitude dispersive 
waves does not apply. A step in this direction is in the work m\mm, where asymptotic stability 
for systems with a single family of nonlinear bound states is treated by a novel time-dependent 
linearization procedure. 


7 A nonlinear toy model of nonlinearity-induced energy transfer 


In this section we explain the idea behind the results on ground state selection and energy equipar¬ 


tition. In the style of sections 5.1 and 5.2 we introduce a toy minimal model, here nonlinear, which 
captures the essential mechanisms and is, by comparison, simpler to analyze. 

Our model is for the interaction of three amplitudes: 


1 . a “ground state” complex amplitude, oto(t) with associated frequency Eq+ 

2 . an “excited state” complex amplitude, aq (t) with associated frequency E\+, and 

3. a “continuum wave-field” complex amplitude, R(x,t), with spectrum of frequencies given by 
the positive real-line [0, oo). 

As with NLS / GP, our model has cubic nonlinearity and the assumption 


cj* = 2£i* — Eq* > 0 


(7.1) 
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is assumed to assure coupling of discrete and continuum modes at second order in the solution 
amplitude. We also introduce a function, x( x )i sufficiently rapidly decaying at spatial infinity. 


(Recall that in section 5.1 we took x( x ) = S(x). 


Our nonlinear model is the following: 


id t a 0 (t) - O 0 q 0 = g (x,R(t)^ aj(t) 
id t ai(t) - fii «i = 2 g (x,R(t)) a 0 (t ) aq(f) 
idt R(x,t) = —A R(x,t) + g x( x ) a i(t) a o(t) 


(7.2) 

(7.3) 

(7.4) 


A first observation is that the system (7.2)-(7.4) is Hamiltonian and has the time-invariant 
quantity 

^ |a 0 (t)| 2 + |«i(t)| 2 + f \R(x,t)\ 2 dx) = 0 


To solve (7.2)-(7.4), we first separate fast and slow scales by introducing slowly varying ampli¬ 


tudes, Aq and A\: 

a 0 (t) = e-^ A 0 (t), a^t) = e " <nit A^t) 

Then, the system for Aq,A\ and R becomes: 

id t A 0 (t) = g (x,R(t)^ A\(t) 

id t Ai{t) = 2 g (x,R(t)) e 1 ^ A 0 (t)Ai(t) 

idt R(x,t) = —A R(x,t) + j xW d((l) 4o(t) 


(7.5) 

(7.6) 

(7.7) 


Note that since the spectrum of —A = [0, oo) and since w* > 0 (by assumption) the Schrodinger 


wave equation, (7.7) is resonantly forced. We proceed now to solve for R(x , t), keeping its dominant 


contributions. Duhamel’s principle and a regularization, motivated by the local decay estimate (7.9) 
below, gives: 


R(t) = e lAt R( 0) - ig e~ I e 

Jo 

pt 

_<5^.0 


AA t 


—i(—A+u*)s 


xA 0 (s) A^s) ds 


= ^u t R _ ige i 


1 


_ d_ 

Jo *(—A + w* — id) ds . 

Next we have, by integration by parts, 


J(-A-cu+)s 


xA 0 (s)A 1 (s) ds 


R(x,t ) ~ e iAt R( 0) - ge~ iu} ** 

_ piAt 

+ j 4o(0)A 2 (0) 

pt 


1 


+ 


—A — w* — iO 

iA(t-s) d 


—A — ca* — 70 

X 


tx X Ao(t)Af(t) 


/ o 


—A — w* — iO ^ ds - 


A 0 (s)Af(s) ds , 
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which may be written as 


m 


-ge 




1 


xx X -Ao{t)Ai(t) + O it 2 


— A — w* — ?'0 

where the error term is obtained using the local energy decay estimate: 


AAt 


(x}~ 


— (x) " = 0(t 2 ) ; t —>• + 00 , 

L 2 (R 3 )-s>L 2 (M 3 ) 


(7.8) 


(7.9) 


—A — w* — iO 

where v is sufficiently large and positive. Substitution of the leading order terms of R(t) in the 


expansion (7.8) into (7.5)-(7.6) and use of the distributional identity 

1 


(-A - w* - iO)" 1 = P.V. 


+ m d(—A — cj*) 


-A - w* 

yields the dissipative system for the amplitudes Ao and A\: 

d t Ao « + r |^| 4 A 0 , dtA! « -2T |A 0 | 2 |Ai| 2 ^i 


(7.10) 


where T ~ g 1 \ x(w*) | 2 > 0 and is generically strictly positive. This is the analogue of the strict 
positivity discussed in the statement of Theorem |6.2| Here, x denotes the Fourier transform of x- 
The omitted correction terms in ( |7.10 ) can be controlled in terms of the quantities ||-R(£)||lp 
(p > 2 and sufficiently large), \\{x)~ v R{t )\\ L 2 and || (x)~ u e zAt (—A — w* — i0) -1 (x)~ u \\ L 2 ^ L 2 . 

Finally, ground state selection with energy equipartition can be seen through the reduction 
(7.10). Indeed, setting Pq = |c^o | 2 , Pi = |au | 2 we obtain 

~ F P 2 P ~ op p 2 p 

~dt rPl Po ’ “ 2rPl Po 

It follows that 2Po(t) + P\(t) ~ 2Po(0) + Pi(0). Sending t — > 00 and using that P\{t) — > 0, (indeed, 
T > 0), we obtain Po(oo) = Pq( 0) + lPi(0). 


8 Concluding remarks 

A fundamental question in the theory of nonlinear waves for non-integrable PDEs is whether and in 
what sense arbitrary finite-energy initial conditions evolve toward the family of available nonlinear 
bound states and radiation. This question is often referred to as the Soliton Resolution Conjecture ; 
see, for example, CEO- The asymptotic stability / nonlinear scattering results of the type discussed 
in section | 6 j being based on normal forms ideas, spectral theory of linearized operators, dispersive 
estimates and low-energy (perturbative) scattering methods, give a local picture of the phase space 
near families of nonlinear bound states. In contrast, integrable nonlinear systems which can be 
mapped to exactly solvable linear motions, allow for a global description of the dynamics mum 
[221:20]. We also note the important line of research by Merle et. al. which is based on the monotone 
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evolution properties of appropriately designed local energies (see, for example, in |79 : , 8 TJ) and does 
not rely on dispersive time-decay estimates of the linearized flow. It would be of great interest 
to adapt these methods for use in tandem with dynamical systems ideas to situations where the 
background linear medium, i.e. a potential V(x), gives rise to structures (defect modes) which 
interact with “free solitons”. 

With a view toward understanding asymptotic resolution in non-integrable Hamiltonian PDEs, 
one of the simplest global questions to ask is the following. Let V denote a smooth, radially- 
symmetric and uni-modal potential well which decays very rapidly at spatial infinity. Assume 
further that — A + V has a finite number of bound states. Consider NLS / GP with a repulsive 
nonlinear potential (g = + 1 ) with radially symmetric initial conditions: 

id t m = - A* + + *(s,0) = *o(M) e H l radial {R d ) 


It is natural to expect that energy which remains spatially localized must reside the family of 
nonlinear defect modes, i.e. bound states, localized within the potential well, which lie on solution 
branches bifurcating from the linear bound states of — A + V [ 88 j; see section [3j Indeed, any 
concentration of energy outside the well should disperse to zero, since for \x\ large the equation 
is a translation invariant NLS equation with repulsive (defocusing) nonlinearity having no bound 
states. We conjecture that for arbitrary initial data solutions either disperse to zero or approach a 
nonlinear defect state, and furthermore that generic initial conditions approach a stable nonlinear 
defect state. 

A step in this direction is the result of Tao 102! for 1 + | < p < 2*, in very high spatial 
dimensions. For spatial dimensions d > 11 there is a compact attractor. That is, there exists K, a 
compact subset of H) adial (R d ), such that K is invariant under the the NLS / GP flow. Moreover, 
i/«e»xR d ) is a global-in-time solution of NLS / GP, then there exists u + E H 1 radial ) 
such that distil i (\l/(i) — e lAt u + ,K ) —> 0, as t —> +oo. 

It is natural to conjecture that K is the set of nonlinear bound (defect) states. 


Finally we mention that the detailed dynamical picture of energy transfer from discrete to 
radiation modes has connections with important questions in applied physics. Note that the Fermi 


golden rule damping matrix, Tq(z,z), appearing in Theorems 6.1 and 6.2, which controls the rate 
of decay of excited states (neutral modes) is controlled by the density of states of the linearized 


operator near the resonant frequency w*; see ( 6 . 2 ). 


Control problem: Can one design the potential, V{x), in NLS / GP in order to inhibit / enhance 
energy transfer and relaxation (decay) to the system's asymptotic state? 

This relates to the problem of controlling the spontaneous emission rate of atoms by modify¬ 
ing the background environment of the atom (for us, the background linear potential), thereby 
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controlling the density of states m A problem of this type was investigated analytically and 
computationally for a closely related parametrically forced linear Schrodinger equation [83]. Such 
a study for NLS / GP could inform the design of experiments, such as those reported on in [78]. 
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